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PRE  FACE 

The  objective  of  the  Atmospheric  Sciences  Laboratory  (ASL)  effort 
within  the  Ground  Target  Signature  (GTS)  program  is  to  develop 
a  total  atmospheric  transmission  model  that  takes  into  account 
molecuiar  and  aerosol  scattering  and  absorption.  The  five- para meter 
generalization  of  the  Zachor  molecular  absorption  model  originally 
proposed  by  Gibson  and  Pierluissi  in  1971  is  in  the  process  of  being 
integrated  into  a  comprehensive  ASL  transmittance  model.  This  absorption 
model  is  reviewed  and  compared  with  the  modified  King's  function 
modei  used  in  the  Air  Force  Cambridge  Research  Laboratory  (AFCRL) 
model  presently  being  used  for  (TTS  purposes.  Since  further  improve¬ 
ments  need  to  be  made  to  develop  the  five-parameter  model  for  general 
use,  additional  work  to  improve  the  model  and  refine  the  computer 
techniques  involved  is  being  done  under  contract  number  DDAD07-73-C- 
0127  at  the  University  of  Texas  at  El  Paso  by  Joseph  Pierluissi, 

Le I  and  Blank,  and  Jerry  Collins.  This  work  i nc I udes  applying  the 
model  to  high  resolution  data,  extending  the  use  of  the  modei  to 
real  atmospheric  conditions,  applying  inhomogeneous  path  techniques, 
and  improving  the  computer  program  efficiency. 

The  authors  gratefully  acknowledge  the  cooperation  and  assistance 
of  Gienn  A.  Gibson  and  Joseph  H.  Pierluissi  of  the  University  of 
Texas  at  Ei  Paso.  We  also  wish  to  acknowledge  the  relevant  discussions 
with  John  E.  Seiby  of  AFCRL  concerning  the  application  of  the  modified 
King's  function  model.  We  especially  wish  to  acknowledge  Richard  B. 

Gomez  for  his  critical  review  and  his  many  helpful  suggestions. 
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INTRODUCTION 


This  is  the  second  in  a  series  of  reports  to  determine  a  total  atmos¬ 
pheric  transmittance  model  for  the  Ground  Target  Signatures  (GTS)  program. 
(Somez  and  Plerluissi  Cl]  reviewed  the  state-of-the-art  of  calculating 
atmospheric  transmission,  and  recommended  that  a  generalized  molecular 
absorption  model  be  considered  as  potentially  the  best  model  available 
for  treating  broad  band  atmospheric  absorption. 

This  report  discusses  and  applies  the  five-parameter  generalization 
of  Zachor’s  model  CZ],  henceforth  referred  to  as  the  FP( f ive-parameter) 
model.  This  application  is  to  a  l-km  horizontal  path  of  water 
vapor  data  calculated  by  Wyatt,  Stull,  and  Plass  C3l  over  spectral 
intervals  of  1200  to  2200  cm~ '  and  4900  to  5800  cm  with  a  resolution 
of  50  cm”  .  To  provide  a  basis  of  comparison,  a  modified  King's 
function  (MKF)  absorption  model,  used  by  Air  Force  Cambridge  Research 
Laboratories  (AFCRL)  in  their  total  atmospheric  transmission  model 
C4],  is  appi ied  to  the  same  path  and  spectral  intervals  using  the 
water  vapor  data  of  Burch.  Singleton,  France,  and  Williams  C5]  degraded 
to  a  resolution  of  50  cm”  . 

The  AFCRL  model  was  chosen  as  a  basis  of  comparison  because  it  is 
used  in  GTS  applications. 

DISCUSSION  OF  THE  F I VE -PARAMETER  MODEL 

Development  of  the  Five-Parameter  Model 

The  FP  band  model  discussed  in  this  paper  is  a  quadratic  generalization 
of  the  four-parameter  absorption  model  developed  by  Zachor  (l2,  6]. 

As  originally  proposed  by  Gibson  and  Pierluissi  [1^,  it  was  developed 
by  considering  the  Zachor  model  as  the  equation  of  an  elliptic 
cone  and  writing  that  expression  as  a  more  general  three-term 
polynomial  equation.  The  Zachor  model  is  a  modification  of  the 
Mayer -(Boody  statistical  band  model  for  gaseous  transmittance. 

The  Mayer-(3oody  model  Is  (I8,  9]: 


T 


exp  — 
.  ( I 


(I) 
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where 


T  =  transmission 

q  =  SU/d  =  kU 

U  =  optical  path  length 

k  =  absorption  coefficient 

S  =  spectral  line  intensity 

d  =  mean  spacing  between  spectral  lines 

3  =  2TrY/d  =  3qP 

3^  =  3  calculated  at  standard  temperature  and  pressure 

Y  =  spectral  line  half-width 

P  =  gaseous  broadening  pressure 

Eq.  (I)  may  be  rewritten  as 


T 


exp 


q3  J 


(2) 


or 


T 


CUP 


(3) 


where  C  =  3Qk  is  the  proportional ity  factor  between  UP  and  q3  and  is 
frequency-dependent.  The  product  CUP  is  used  instead  of  q3  because  q 
and  3  cannot  be  directly  measured. 
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For  U/P<<l,  the  weak-line  approximation,  Eq .  (3)  reduces  to 


T  =  exp(-kU) 
w  ^ 


(.4) 


and  forU/P»l,  the  strong-line  approximation,  Eq .  (3)  reduces  to 


Ts  =  exp 


CUP.Js  j 


-  (^) 


(5) 


Combining  Eqs.  (3),  (4),  and  (5),  the  Mayer-Goody  model  may  be  re¬ 
written  as 


9  -2  “2 

( lnT)“  =  ( InT  )  +  (InT  ) 

W  s 


C6) 


Zachor  C2l  noted  that  this  equation  may  be  regarded  geometrically 
as  the  resultant  magnitude  of  two  perpendicular  vectors  of  lengths 
(In  T  )“  and  (In  Tj)"',  es  shown  in  Figure  1(a).  Extending  this 
geomeTric  interpretation  and  applying  the  cosine  law,  he  proposed 
a  more  general  model  which  may  be  expressed  as 


Z2  =  (^)  x2  +  y2  -  (M)  xy 


(7) 


where 

Z  =  I/In  T 

X  =  l/U 

y  =  I/I n  T 

s 

M  =2  cos  6 
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(  JUvTs)-l 


Figure  la.  The  vector  representation  of  the  Mayer-Goody  model. 


Figure  lb.  The  vector  representation  of  the  Zachor  model. 


R 


(In  t)  Is  now  the  sum  of  two  vectors,  (In  t^)  and  (In  t  )  ,  whose 

directions  differ  by  a  fixed  angle  6,  as  shown  In  Figure  1(B).  The 
value  6  has  no  effect  on  the  value  of  t;  however,  it  does  control  the 
rate  of  transition  between  t  and  .  The  cross  term  can  thus  be 
Interpreted  as  representing  The  region  of  intermediate  absorption. 

The  Zachor  formula,  when  written  in  the  form  of  Eq.  (7),  is  the 
equation  of  an  elliptic  cone.  The  ellipse  is  in  the  xy-plane,  and 
the  xy  term  corresponds  to  a  rotation  of  the  ellipse  about  the  z-axis. 
Since  M  =  2  Oos  6,  M  will  be  less  than  or  equal  to  2.  This  restriction 
guarantees  that  the  cone  will  be  elliptic  (or  parabolic  for  |m|  =  2). 

This  comes  from  the  restriction  In  analytic  geometry  that  <  4A|A2, 

where  A|  is  the  coefficient  of  x^,  A2  Is  the  coefficient  of  y^,  and 
Aj  is  the  coefficient  of  xy. 

To  use  the  analytical  models  described  above,  one  essentially  de¬ 
termines  an  elliptic  cone  that  best  fits  the  data  points  (xj,  yj,  Zj). 
This  may  be  done  by  using  minimization  techniques  to  weight  the  terms 
of  the  equation  of  the  ellipse.  The  coefficients  are  determined  from 
the  resulting  best  fit.  This  procedure  must  then  be  repeated  for 
each  frequency,  since  the  parameters  are  frequency-dependent. 

Using  the  observations  presented  above,  Gibson  and  Fieri uissi  C7] 
proposed  a  generalized  polynomial  model  which  would  allow  for  complete 
flexibility  in  weighting  the  terms  of  the  equation.  The  proposed  model 
i  s 


Z2  =  B  X2  +  B  +  B  XY 
w  s  ws 


(8) 


where  B^  is 
tran 


the  frequency-dependent  parameter  which  weights  the  weak- 
mission  with  respect  to  the  other  two  terms.  It  is  not  necessar- 

is  the  frequency-dependent 


line 

ily  (l/K^),  as  in  the  Zachor  model .  B^ 
parameter  which  weights  the  strong-line  transmission  terms.  B^g  is 
the  frequency-dependent  parameter  which  weights  the  transition  between 
X  and  y.  These  coefficients  satisfy  the  condition  B^g<4B^yB5. 


In  a  later  paper  Pierluissi  ClOl]  further  generalizes  the  model  and 
describes  the  general  polynomial  transmission  model  as 
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+  . .  . 


C9) 


Z2  =  B  X2  +  B  y2  +  B^^XY  +  B,  X^Y  +  B^^XY^ 
w  s  ws  3ws  2WS 


in  which  the  right  side  of  the  equation  is  simply  an  n-order  poly¬ 
nomial  given  by  the  expression 


F(x,  y) 


L. 


B  x'^y' 

n . 

I 


(  10) 


The  first  order  terms  of  Eq .  (10)  correspond  to  a  trans I  at i on >of  the  x- 
and  y-axes  and  may  be  removed  by  algebraic  methods.  The  accuracy 
produced  by  the  third  order  and  higher  terms  is  not  needed  for 
general  atmospheric  transmittance  modeling,  and  the  inclusion  of  these 
terms  results  in  greater  complexity  of  the  model  and  more  computational 
effort.  Eq .  (9)  is  therefore  truncated  after  the  first  three  terms, 
yielding  Eq .  (8).  As  more  terms  are  used,  the  coefficients  lose 
their  physical  meaning. 

As  written,  the  Zachor  model,  Eq .  (7),  and  the  FP  model,  Eq .  (8), 
have  only  two  and  three  explicit  parameters,  respectively.  However, 
both  models  use  the  two-parameter  strong- line  transmission  function 
proposed  by  King  CM]  for  their  representations  of  t  .  King's 
formu la  Is 


T 

S 


I  -  P{n,  Cnr(n) 


(II) 


where  P(a,  x)  is  the  incomplete  gamma  function 


10 


(.12) 


P(a,  x)  =  [r(a)I| 


-i 


/ 


,3-  I 


dt 


King  developed  his  strong-line  transmission  model  based  on  the  fact 
that  the  character  of  the  line  spacing  is  the  most  important  factor 
in  the  transmittance  of  overlapping  spectral  lines.  He  compared 
the  transmittance  of  the  Elsasser  model  to  the  transmittance  of  the 
Mayer-(3oody  model  and  found  no  relationship  between  the  two;  however, 
in  taking  the  derivatives  of  the  absorptances  of  the  models  with 
respect  to  the  absorber  quantity,  he  discovered  a  relationship 
which  led  him  to  formuiate  two  important  assumptions;  (I)  The 
absorption  derivative  of  widely  spaced  overlapping  Lorentz  lines 
can  be  expressed  as  the  product  of  the  absorption  derivative  of 
the  lines  considered  as  non-interacting  absorbers  and  an  overlap 
factor  which  involves  the  interaction  of  the  neighboring  lines.  (2) 

For  widely  spaced  Lorentz  absorption  lines,  the  overlap  factor  can  be 
expressed  as  an  exponential  function  In  which  the  power  of  the  argument 
is  related  to  the  variance  of  the  line  spacing.  The  assumptions 
form  the  basis  for  the  derivation  of  King's  formula. 

Eg.  (II)  has  two  frequency-dependent  parameters,  n  and  C.  The 
parameter  n  is  an  adjustable  strong-line  parameter  which  controls  the 
ratio  of  the  variance,  ,  and  the  squared  mean  of  the  line  spacing, 
d^,  for  widely  spaced,  equally  Intense  strong  lines,  as  in  the 
probability  density  function  L6!]: 


a 


2 


d^ 


{2nr(2n)/[nr(n)]2}-| 


(  13) 


The  Mayer-(3oody  model  (Poisson-distributed  lines)  is  characterized  by 
unit  variance,  while  the  Elsasser  model  (regularly  spaced  lines)  is 
characterized  by  zero  variance.  For  Poisson-distributed  lines  (n  =  I), 
Eq.  (13)  gives  the  correct  result;  however,  for  regularly  spaced  lines 


I  I 


(n  =  ^) ,  It  gives  a^/d^  =  (4/ir-l)  Instead  of  zero.  Nevertheless,  King’s 
formula,  Eq .  (II),  is  very  useful  because  it  provides  a  continuous 
set  of  trial  functions  for  fitting  strong-line  data,  and  it  reduces  to 
the  Elsasser  strong-line  approximation  for  n  =  h  and  to  the  Mayer-(3oody 
strong  line  approximation  for  n  =  I.  For  n  >  I,  the  expression 
represents  clustering  of  spectral  lines.  Eq .  (II)  is  thus  a  generalized 
strong-line  absorption  model  whose  line  spacing  goes  from  complete 
regularity  to  complete  randomness  and  to  the  limit  of  clustering  in 
which  the  spectral  lines  are  superimposed  upon  one  another. 

The  parameter  C  is  a  strong- line  parameter  related  to  the  line 
spacing,  d,  the  line  strength,  S,  and  the  half-width,  y  ,  of  the 
spectral  lines  by 


2tty  S 

C  =  - ^  (  14) 

d^ 


With  King's  stnong-line  transmission  function  (Eq.  (II))  incorporated 
into  Eq.  (8),  the  proposed  model  becomes  a  five-parameter  model 
(n,  C,  Bg,  All  five  parameters  must  be  calculated  for  each 

frequency.  The  techniques  used  to  determine  these  parameters  will 
be  discussed  in  the  next  section. 

That  the  FP  model  is  indeed  a  very  general  and  very  useful  analytical 
model  can  be  shown  in  that  at  least  nine  other  models  can  be  derived 
from  it  by  forcing  the  parameter  to  certain  values  Cl2ll. 

The  models  that  can  be  obtained  are: 

2 

1.  Beer's  Law.  By  setting  B  =  B  =0  and  B  =  l/k  ,  we  obtain 

x  =  -  ^  s  ws  w 

2.  King's  Strong-1 i ne  Model .  By  setting  B  =  B  =0  and  B  =  I, 
we  obtain  Eq.  (II),  repeated  here  for  convenience. 


T  =  l-P{n,rCnT(n) 

S  IT 
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3.  Elsasser’s  Strong-line  Approximation.  By  setting  B  =  B  =0, 

B  =1,  and  n  =  0.5,  we  obtain 
s  ' 


T 


vlh,  ( 


CUP^^ 

2 


] 


=  I  -  erf 


(CUP)J5 


(  15) 


4.  The  Mayer-Goody  Strong- line  Approximation.  By  setting  B  = 

B  =  0,  B  =  I and  n’  '=  I,  we  obtain 
ws  s 


p(i,  (22^)'^) 


=  exp  C-( 


tr 


5.  The  Mayer-Goody  Model .  By  setting  B^^  =  0,  B^ 
and  n  =  I  ,  we  obtain  Eg.  (3),  repeated  here. 


(  16) 


exp  C-  ( 


rkUF 


CUP 


6.  The  Generalized  Mayer-Goody.  By  setting  B  =0  and  using 
King’s  expression  for  with  n  variable,  we  obtain 


( Inx ) 


I 

(kUP 


+  (Inx  ) 
s 


(  17) 


7.  The  Modified  Elsasser.  Setting  B  =0  and  n  =  *5,  we 
obtain 
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(kU)^ 


(  18) 


( Int) 


-2 


+  Cl-erf( 


CUP,isn-2 
2  ^ 


„  8.  Zachor’s  Model  .  By  setting  B 
l/k  ,  we  obta  i  n  ^ 


I ,  B  =  M/k,  and  B  = 
'  ws  '  w 


(  Int) 


-2 


I 

TWP 


+ 


(Inis 


-2 


_M 

kU 


(  I nx 

s 


-I 


(  19) 


9.  Intermediate  Absorption  (any  model).  Setting  B  =  B  =  0, 
we  obtai n 


T  =  B,,^lnT  Inx 
I  ws  w  s 


(20) 


Since  these  other  models  can  be  derived  from  the  FP  model,  once  the 
FP  model  is  computerized,  one  can,  with  appropriate  fixing  of  parameters, 
obtain  any  of  the  above  models  for  use  in  the  calculation  of  gaseous 
transmittance.  A  reduction  in  the  number  of  parameters  will  result 
In  a  decrease  of  computer  time  and  of  required  oomputer  storage;  however. 
It  will  also  result  in  a  decrease  in  the  accuracy  of  the  model. 
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Mathematical  Techniques  Used  to  Determine  the  Five  Parameters 

The  mathematical  algorithms  used  by  Gibson  and  Pierluissi  to  determine 
the  five  parameters  will  be  discussed  in  this  section.  As  mentioned 
in  the  previous  section,  the  parameters  are  frequency-dependent  and 
must  be  calculated  for  each  frequency.  Ideally,  Eq .  (8)  should  be 
minimized,  solving  for  all  five  parameters  simultaneously  by  using 
al i  the  absorption  data;  however,  such  a  technique  has  not  been 
perfected. 

A  procedure  similar  to  that  used  by  Zachor  C2]  was  used  to  determine 
the  parameters.  First  the  strong- line  parameters,  n  and  C,  were 
determined.  Then  these  values  were  used  in  calculating  the  quadratic 
parameters:  B  ,  B^,  and  B  .  A  simplified  flowchart  illustrating  the 
computer  procedures  described  in  this  section  is  shown  in  Figure  2. 

The  procedure  for  determining  n  and  C  is  based  on  the  fact  that  for 
a  fixed  UP  product,  such  as  U/P>>l  (the  strong-line  approximation),  the 
absorptance  must  approach  an  upper  limit  of  I-T3,  where  T3  Is 
the  strong-line  transmittance.  The  strong-line  parameters  can  thus 
be  determined  from  the  experimental  data  by  plotting  log  (I-T3)  versus 
log  (UP)  and  selecting  from  these  data,  points  which  are  dominated  by 
strong- line  absorptance.  These  points  are  then  fitted  to  a  surface 
obtained  by  plotting  the  log  of  the  theoretical  strong-line 
transmission  function,  T3,  given  by  Eq .  (II),  versus  the  log  of  the 
quantity  (CUP)  versus  the  parameter  n.  The  value  of  Log  C  may  be 
obtained  by  taking  the  difference  between  Log  (CUP)  and  Log  (UP)  from 
these  two  piots.  The  value  of  n  may  be  obtained  directly  from  the 
displacement  along  the  n-axis.  The  correct  values  of  n  and  C  are 
obtained  when  the  following  expression  is  minimized: 


i=N 


(21) 
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Figure  2.  Flowchart  of  the  five-parameter  model. 
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Figure  2  (con . ) 
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where  D=Log  C,  Pj  =  log  (  I -x  | ) ,  f(n,  D)  =  log(P{n,Cnr(n)  (2CUP/tt)^I]'^^), 
and  N  =  the  number  of  data  points. 

The  data  used  in  the  calculation  of  the  strong-1 ine  parameters  were 
computer-selected  by  setting  four  levels  of  absorptance:  L| ,  L2, 

L-r,  and  14*  Levels  L|  and  L2  were  used  to  determine  which  data  would 
be  used  in  the  strong-line  calculations  (see  Figure  2).  Levels  L^  and 
L4  were  set  to  exclude  values  of  absorptance  near  0^  and  100^, 
which  in  general  are  questionable.  Small  differences  in  these  values 
may  produce  large  errors  in  the  parameters  to  be  calculated.  These 
levels  were  predetermined  by  trial  and  error,  and  it  was  found  that 
the  levels  L|  =  0.5,  L2  =  0.2,  L-^  =  0.005,  and  L4  =  0.995  produced 
the  best  parameters  for  water  vapor. 

First  the  data  were  checked  against  L.;  no  values  greater  than  L^ 

were  used.  The  data  were  next  checked  against  L|,  and  if  there  were  at 

least  three  data  points  greater  than  L|,  then  all  of  those  points 

were  used  in  the  calculations.  If  there  were  not  three  points  above 

L|,  then  the  data  were  checked  with  L2*  If  there  were  three  or  more 

points  above  L2,  then  those  points  were  used;  n  was  then  set  equal 

to  a  predetermined  set  value  and  Eq .  (14)  was  minimized  with  respect 

to  D  only.  This  was  necessary  since  minimization  with  respect  to 

n  and  C  was  found  to  be  difficult  under  these  circumstances  Cl33. 

The  selection  of  a  value  of  n  was  made  by  a  trial  and  error  procedure, 

and  a  value  of  one  was  found  to  produce  good  parameters.  If  there  were 

only  one  or  two  points  above  L2,  then  the  three  highest  absorptance 

values  were  used  for  the  minimization  with  respect  to  D  only.  If 

there  were  no  points  above  L2,  then  the  three  highest  absorptance 

values  were  used  for  the  minimization  with  respect  to  D  only.  If 

there  were  no  points  above  L2,  the  absorption  was  assumed  to  be 

weak- I ine.  The  points  between  L3  and  L2  were  then  selected  and  Beer's 

Law  was  used.  B  and  B  were  then  set  to  zero  and  only  B  was 
,  ,  ,  .  s  ws  '  w  ' 

ca I cu I ated . 

Since  Tg  Is  nonlinear  in  n  and  D,  the  procedure  used  to  minimize 
Eq.  (14)  was  an  iterative  nonlinear  technique.  A  modified  Newton- 
Raphson  method  was  used  Cl 43;  this  expression  was  of  the  form 


Vl 


g<=k> 


(22) 
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where  g  (aj^)  Is  the  set  of  equations  obtained  by  taking  the  partial 
derivatives  of  Eq.  (14)  with  respect  to  n  and  C,  respectively,  and 
setting  them  equal  to  zero,  and  is  the  set  of  equations  obtained 
by  taking 


8n 


(23) 


and 


8g(a^) 

Id 


(24) 


The  a, 's  are 
k 


a 


k 


k  =  0, 


I, 


(25) 


The  matrix  obtained  for  g(a|^)  is 


= 


i  =  N 
r' 
4. 

i=l 


9f  , 

9D  ^i 


f ) 


i=N 

I 

i=l 


i 

ii 

9n 


(p  j 


f) 


k  =  0,  I ,  . . . 


(26) 
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The  matrix  M|^  is  obtained  by  expanding  the  function  f(n,  D.)  in 
the  Taylor  series  about  the  point  (n^  Dj),  which  makes  f  a  minimum. 

The  function  can  be  closely  approximated  in  the  region  of  interest 
by  considering  only  the  linear  terms  in  the  Taylor  series  expansion  of 
f.  Using  this  simplifying  approximat ion^  we  obtain 


(27) 


where  the  symbol  indicates  that  the  expression  is  evaluated  at 
the  point  (n.,  ^j)*  The  technique  involved  requires  making  an  initial 
estimate  of  the  values  of  n  and  D,  then  making  successive  guesses 
by  incrementing  n  by  no  more  than  0.1  and  D  by  no  more  than  0.3 
until  successive  are  found  such  that  their  elements  are  less 

than  0.0005.  When  this  difference  is  obtained,  Eq .  (14)  is  minimized 
with  respect  to  n  and  D.  In  this  procedure  n  is  minimized  first, 
and  this  value  is  used  in  minimizing  D.  In  order  for  n  and  D  to 
be  minimized  in  this  fashion,  the  initial  estimates  must  be  such 
that  n  and  D  can  be  determined  within  a  limited  number  of  iterations, 
which  was  arbitrarily  chosen  as  46  in  the  five-parameter  computer 
program. 
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For  the  cases  when  n  was  set  equal  to  1.0  and  Eq .  (14)  was  minimized 
with  respect  to  D  only,  Eqs.  (26)  and  (27)  became 


i=N  i 


3‘°k>  =  Z  I?  ‘"i''’ 


k  =  0,  I 


(28) 


and 


i  =  N 


M 


k 


z 


ilf 

ao  / 


k  =  0,  I 


(29) 


In  order  to  determine  the  matrices  for  M|^  and  g(a|^)^  +he  function 
f(n,  D)  and  its  partial  derivatives  with  respect  to  n  and  D  must 
be  taken  and  evaluated.  A  series  approximation  method  may  be 
used  to  evaluate  the  incomplete  gamma  function;  however,  the  function 
converges  very  slowly  when  it  is  close  to  the  value  I.  For  this 
reason  the  (n,  UP)  plane  was  divided  up  into  several  areas  of 
50  points  by  50  points.  To  approximate  the  log  of  the  incomplete 
gamma  function,  a  set  of  general  third-order  polynomials  In  n  and 
log  (CUP)  was  written  in  the  form 


f(n,  x)  =  A.x^  +  A_nx^  +  A,n^x  +  A.n^  +  A 
12  3  4  5 


(30) 


where  x  =  log  (CUP).  Tabulated  values  for  the  Incomplete  gamma  function 
were  used  Cl5ll,  and  Eq .  (30)  was  fitted  to  those  values  by  using 
least-squares  techniques  to  determine  the  coefficients.  These 
polynomials  and  their  part ia Is  with  respect  to  n  and  D  were  used 
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in  the  calculation  of  g(aj^)  and  A  similar  polynomial  expression 

for  the  incomplete  gamma  function  was  also  determined,  since  it 
must  be  evaluated  in  order  to  find  values  for  The  functions 

were  evaluated  in  the  regions  0.5^n^3.0  and  -4^log  (CUP)^4. 

The  polynomial  approximations  were  within  0.001  of  the  actual  values 
of  the  functions. 

After  the  parameters  n  and  C  were  determined,  the  quadratic  parameters 
0^,  Bg,  and  were  obtained  by  applying  least-squares  techniques 
to  minimize  the  weighted  difference  between  the  theoretical  and  the 
experimental  i.e.. 


i  =  N 

'y  w(t^)CB  x.^  +  B  y.^  +  B  x.y.  - 

^  •  I  W  I  S  I  ws  I  '  I  I 

i  =  l 


(31) 


2  6 

where  wCtp  =  (  Intp  ,  and  it  is  the  measured  value  of  the  trans¬ 

mission.  This  weighting  function  was  postulated  from  the  assumption 
that  the  transmittance  is  a  random  variable  with  a  probability  density 


S  =  ( Ini ) 


(32) 


2  6  —  I 

The  variance  of  this  density  is  proportional  to  Ct*  (Ini*)  ]**  . 

In  least-squares  methods  the  weighting  function  is  conventionally 
inversely  proportional  to  the  variance  CIS].  The  expression  for 
w(t*)  Is  thus  Justified  on  this  premise. 

The  procedure  for  minimizing  Eq .  (31)  involves  taking  the  partial 
derivatives  with  respect  to  B^,  Bg,  and  B^g  and  setting  them  equal 
to  zero.  This  yields  three  equations  to  be  solved  for  these 
parameters : 


i=N 

'y  w(t^)X.2 (X.2b  + 
Z  .  II  I  w 


^i^iEws  +  Y.2B 

I  I  IS 


-  Z*2)  =  0 

I 


(33a) 


(33b) 


>: 

i=i 


w(t*)Y.2(X.2b  + 

II  I  W 


X-  +  Y.2B 

I  I  ws  IS 


_  Ztf2)  =  0 
I 


i=N 

/  w(t^)X.Y.  (X.2b  +  X.Y,B^-  +  Y.2b  -  Zf2)  =  0  (33c) 

/  ,  I  I  I  I  W  I  i  ws  IS  I 


The  value  of  Xj  Is  determined  from  l/U  and  the  value  of  yj  is 
determined  from  the  third-degree  polynomial  approximation  for  the 
Incomplete  gamma  function.  The  above  set  of  equations  Is  easily 
solved  at  each  frequency  for  the  three  quadratic  parameters. 

APPLICATION  TO  WATER  VAPOR  DATA  AND  COMPARISON  TO  AFCRL  MODEL 

The  FP  model  was  applied  to  the  water  vapor  absorption  data  of  Wyatt, 
Stull,  and  Plass  [13]  over  absorption  bands  of  1 200  to  2200  cm“ '  and 
4900  to  5800  cm“ '  with  a  resolution  of  50  cm  .  The  transmi ttances 
for  the  l-km  horizontal  atmospheric  path  were  calculated.  The 
selection  of  this  path  resulted  in  an  equivalent  water  vapor 
concentration  of  2  pr-cm.  The  results  of  this  application  are 
represented  by  the  solid  lines  In  Figures  3  and  4. 

The  AFCRL  model,  described  by  McClatchey  et  al  .  C4l],  employs  a 
modified  King's  function  (MKF)  as  its  molecular  absorption  model. 

It  is  an  empirical  model  specifically  tailored  to  fit  a  specific  set 
of  data  for  water  vapor,  carbon  dioxide,  ozone,  and  a  combination 
of  uniformly  mixed  gases  (oxygen,  methane,  carbon  monoxide,  and 
n itrogenoxi de) .  Since  the  AFCRL  model  has  found  use  in  GTS  applications, 
a  comparison  between  It  and  the  FP  model  was  made. 

King  dll]  attempted  to  write  a  general  expression  for  atmospheric 
transmittance  (t)  which,  under  certain  conditions,  would  approach 
either  the  strong-line  or  the  weak-line  approximation  of  either  the 
Elsasser  model  or  the  Mayer-(3oody  model.  He  proposed  that  t  take 
the  general  functional  form 
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Figure  3.  Comparison  of  the  FP  model  with  the  AFCRL  absorption  model 
over  the  spectral  interval  1200-2200  cm"'  with  a  resolution  of  50  cm" 


Figure  4.  Comparison  of  the  FP  model  with,  the  AFCRL  absorption  model 
over  the  spectral  interval  4900-5800  cm"  with  a  resoiution  of  50  cm" 
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(34) 


T  (u)  =  f[C(u)  alp"] 
Av 


where 

T^^(u)  =  transmittance  averaged  over  the  spectral  Interval 

C(u)  =  frequency-dependent  absorption  coefficient 

AL  =  optical  path  length 

P  =  effective  broadening  pressure 

n  =  frequency-dependent  parameter 


The  King  model  Is  therefore  a  two-parameter  model,  with  n  and  C 
as  the  frequency-dependent  parameters.  In  the  case  of  the  MKF 
model,  a  mean  value  for  n  was  determined  graphically  for  a  range  of 
frequencies,  thus  reducing  It  to  a  one-parameter  model  Cl 7]. 


In  order  to  compare  equivalent  water  concentrations  for  a  I -km 
horizontal  path,  the  AFCRL  tropical  atmosphere  model  was  chosen. 

The  results  are  represented  by  the  dashed  lines  In  Figures  3  and  4. 


In  the  1200-2200  cm  band  In  Figure  3,  we  note  that  the  two  models 
predict  generally  similar  though  not  Identical  results.  In  the 
region  between  1850-2200  cm”  ,  the  two  models  differ  significantly. 
Figure  4  the  two  models  differ  In  the  region  4900-5100  cm” L  In 


In 

order  to  understand  this  difference,  the  Input  data  must  be  examined 


Synthesized  Input  data  and  experimental ly  measured  data,  both  from 
Burch  et  al .  CS],  were  used  In  the  development  of  the  MKF  model. 

The  da+a  of  i^yatt  et  a  I  .  C3]  were  obtained  theoretical  ly  by  calculating 
the  absorptions  due  to  spectral  lines  corresponding  to  transitions 
between  various  vibration-rotation  energy  levels.  Discrepancies 
between  experimental  and  theoretical  data  may  be  due  to  Wyatt's 
application  of  the  equivalent  symmetric  rotor  approximation  to  the 
highly  asymmetric  water  molecule  C4],  since  in  the  theoretical  cal¬ 
culations  it  was  necessary  to  assume  that  the  water  molecule  was  a 
symmetric  rotator.  Both  sets  of  data  are  shown  in  Figures  5  and_6 
(obtained  from  Wyatt  et  a  I . )  for  1200-2200  cm”  and  4900-5800  cm  , 
respective! y . 
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Figure  5.  Comparison  of  the  theoretical  calculations  of  the  transmittance 
of  the  6.3  y  band  with  the  experimental  measurements  of  Burch  et  al. 

The  theoretical  values  have  been  averaged  over  a  20  cm“ '  i nterva I .* 

*  Originally  Figure  I  from  C3l]. 
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Figure  6.  Comparison  of  the  theoretical  calculations  of  the  transmittance 
of  the  1.87  y  band  with  the  experimental  measurements  of  Burch  et  a  I . 

The  theoretical  values  have  been  averaged  over  a  20  cm“ ’  interval.t 
t  Originally  Figure  2  from  [33. 
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Figure  7.  Comparison  of  the  FP  model  to  theoretical  data  of  Wyatt  et  al . , 
and  comparison  of  the  AFCRL  rrode  I  to  empirical  data  of  Burch  over  the 
spectral  interval  1200-2200  cm"'  with  a  resolution  of  50  cm"'. 
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Figure  8.  Comparison  of  the  FP  model  to  the  theoretical  data  of  Wyatt 
et  al.,  and  a  comparison  of  the  AFCRL  model'  to  the  empirical  data  of 
Burch  over  the  spectral  interval  4900-5800  cm"'  with  a  resolution  of 
50  cm" ' . 
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A  second  comparison  was  then  made  between  the  results  of  the  models 
and  their  respective  data.  The  conditions  for  this  comparison 
are  as  shown  in  Figures  5  and  6.  The  results  are  shown  in  Figures 
7  and  8,  where  the  solid  line  represents  the  results  of  the  FP  model 
and  the  dashed  line  that  of  the  MKF  model.  The  dots  show  the 
original  Wyatt  et  al .  data  points,  and  the  squares  the  published 
data  of  Burch  aj_.  Tables  I  and  2  show  the  absolute  deviation 
and  give  the  RMS  deviation  of  all  points  over  each  interval  for 
both  models.  The  RMS  deviation  for  the  1200-2200  cm“ '  band  is  1.9 
X  10“^  for  the  FP  model  and  71 .3  x  10“^  for  the  MKF  model;  for  the 
4900-5800  cm“ '  band,  it  is  I .7  x  10“^  for  the  five-parameter  model 
and  48.2  x  10“^  for  the  MKF  model.  The  FP  model  reproduced  its 
data  within  0.5^  accuracy,  while  the  MKF  model,  with  the  exception 
of  three  points,  reproduced  the  pub  I  ished  data  of  Burch  et  a  I . 
to  within  10/^. 


CONCLUDING  REMARKS 

The  addition  of  a  fifth  parameter  to  generalize  the  Zachor  model, 
combined  with  the  employment  of  high  speed  digital  computers  to 
accurately  determine  the  parameters  of  the  model  by  least-squares 
techniques  and  non-linear  iteration  methods,  makes  the  five-parameter 
molecular  absorption  model  a  highly  accurate  and  the  most  general 
analytical  formula  currently  available.  In  contrast  to  empirical 
models,  the  f i ve- parameter  model  also  has  the  advantage  of  not 
being  restricted  to  any  particular  set  of  data.  Consequently, 
the  five-parameter  molecular  absorption  model  is  to  be  integrated 
into  the  Atmospheric  Sciences  Laboratory's  (ASL)  total  atmospheric 
transmittance  model . 

In  the  1200-2200  cm  '  band,  the  five-parameter  model  and  the  modified 
King's  function  model  produce  comparable  though  not  identical 
results;  in  the  1850-2200  cm“ '  region,  they  differ  markedly.  In 
the  4900-5800  cm“ '  band  they  also  predict  similar  results  except  in 
the  4900-5100  cm~ '  region,  where  they  again  differ  markedly.  The 
results  of  the  comparison  between  the  five-parameter  model  and  its 
input  data  and  the  modified  King's  function  model  and  the  published 
data  of  Burch  et  al .  indicate  that  the  former  more  accurately 
reproduces  its  data.  To  provide  a  more  specific  comparison,  the 
two  models  need  to  be  developed  and  appi ied  to  the  same  set  of 
data  over  various  atmospheric  paths. 
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TABLE  I 


CX3MPARIS0N  OF  THE  PREDICTIONS  OF  THE  FIVE-PARAMETER  ABSORPTION  MODEL 
AND  THE  MODIFIED  KING'S  FUNCTION  ABSORPTION  MODEL  WITH  THEIR  RESPECTIVE 
INPUT  DATA  OVER  THE  1200-2200  cm"'  WATER  VAPOR  BAND 


Wave 

Number 

Wyatt 
et  al  . 
Data 

FP 

Model 

Absol ute 

Dev  i at  ion 

Burch 

Data* 

MKF 

Mcdel * 

Absol ute 
Deviat ion* 

1200 

.972 

.975 

.003 

.98 

.96 

.02 

1250 

.954 

.955 

.001 

.94 

.91 

.03 

1300 

.781 

.779 

.002 

.87 

.83 

.04 

1350 

.518 

.521 

.003 

.60 

.50 

.  1  0 

1400 

.299 

.302 

.003 

.06 

.20 

.14 

1450 

.  1  55 

.  1  55 

.000 

.10 

.13 

.03 

1500 

.047 

.046 

.001 

.01 

.01 

.00 

1  550 

.118 

.114 

.004 

.00 

.01 

.01 

1600 

.165 

.  165 

.000 

.12 

.36 

.24 

1650 

.073 

.070 

.003 

.00 

.01 

.01 

1700 

.052 

.051 

.001 

.01 

.01 

.00 

1750 

.082 

.083 

.001 

.02 

.03 

.01 

1800 

.  123 

.  122 

.001 

.10 

.16 

.06 

1850 

.289 

.290 

.001 

.30 

.37 

.07 

1  900 

.435 

.437 

.002 

.52 

.58 

.06 

1950 

.557 

.559 

.002 

.65 

.64 

.01 

2000 

.691 

.692 

.001 

.77 

.77 

.00 

2050 

.773 

.773 

.000 

.89 

.88 

.01 

2100 

.850 

.849 

.001 

.94 

.94 

.00 

21  50 

.915 

.914 

.001 

.96 

.94 

.02 

2200 

.957 

.958 

.001 

.99 

.98 

.01 

FP  Model  RMS  Deviation  =  1.9  x  10 


MKF  Model  RMS  Deviation  =  71.3  x  10 

*On I y  two  significant  figure  accuracy  was  available  for  these  data. 


29 


TABLE  2 


COMPARISON  OF  THE  PREDICTIONS  OF  THE  FIVE-PARAMETER  ABSORPTION  MODEL 
AND  THE  MODIFIED  KING'S  FUNCTION  ABSORPTION  MODEL  WITH  THEIR  RESPECTIVE 
INPUT  DATA  OVER  THE  4900-5800  cm“ '  WATER  VAPOR  BAND 


Wave 

Number 

Wyatt 
et  a  1  . 
Data 

FP 

Model 

Abso 1 ute 

Dev i at  ion 

Burch 

Data* 

MKF 

Model  * 

Absol ute 
Devlat ion* 

4900 

.978 

.979 

.001 

1  .00 

1  .00 

.00 

4950 

.976 

.976 

.000 

1  .00 

.99 

.01 

5000 

.965 

.965 

.000 

.99 

.97 

.02 

5050 

.934 

.935 

.001 

.97 

.96 

.01 

5100 

.858 

.856 

.001 

.90 

.85 

.05 

5150 

.720 

.717 

.003 

.71 

.67 

.04 

5200 

.571 

.572 

.001 

.59 

.45 

.04 

5250 

.554 

.556 

.002 

.42 

.38 

.04 

5300 

.41  1 

.407 

.004 

.30 

.22 

.08 

5400 

.520 

.521 

.001 

.44 

.41 

.03 

5450 

.599 

.599 

.000 

.40 

.35 

.05 

5500 

.679 

.680 

.001 

.63 

.51 

.14 

5550 

.770 

.777 

.002 

.81 

.76 

.05 

5600 

.899 

.899 

.000 

.94 

.91 

.03 

5650 

.970 

.971 

.001 

.97 

.97 

.00 

5700 

.993 

.994 

.001 

.99 

.97 

.02 

5750 

.995 

.997 

.002 

1  .00 

.98 

.02 

5800 

.998 

.998 

.000 

1  .00 

.99 

.01 

FP  Model  RMS  deviation  =  I .7  x  10 


MKF  Model  RMS  deviation  =  48.2  x  10 

^On I y  two  significant  figure  accuracy  was  available  for  these  data. 
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Currently  developmental  work  is  being  done  for  the  ASL  GTS  modeling 
program  (Contract  Number  DDAD07-73-C-0I 27)  at  the  University  of 
Texas  at  El  Paso  to  extend  the  five-parameter  model  to  real 
atmospheric  conditions.  This  work  consists  of  the  application 
of  inhomogeneous  path  techniques  to  the  transmission  calculations 
and  the  inclusion  of  transmission  parameters  for  all  important 
gaseous  atmospheric  constituents.  The  computer  program  efficiency 
will  be  maximized  and  the  model  will  be  applied  to  high  resolution 
data. 

Because  of  the  large  number  of  computations  necessary  to  compute  the 
five  parameters,  the  model  is  best  used  over  a  limited  spectral 
interval.  However,  if  high  accuracy  is  not  required  or  if  computer 
time  and  storage  is  limited,  the  completely  developed  program 
will  have  the  capability  of  setting  one  or  more  of  the  parameters 
to  a  fixed  value.  This  will  result  in  reducing  the  five-parameter 
model  to  one  of  the  nine  models  mentioned. 
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Dr.  A.  D.  Belmont 
Research  Div 
Control  Data  Corp 
Minneapolis,  Minn  65440 

Chief,  Aerospace  Environ  Div 
CODE  S&E-AERO  Y 
NASA 

Marshall  Space  Flight  Center,  Alabama  36802 
Commandant 

US  Army  Field  Artillery  School 
ATTN:  ATSFA-G-RA.RK 
Ft  SiU,  OK  73503 

Technical  Library 
WSMR,  New  Mexico  88002 
3 


Natl  Aeronautics  &  Space  Admin 
George  C  Marshall  Space  Fli^t  Center 
Aero-Astrodynamics  Laboratory 
ATTN:  S&E-AERaVS/Mr.  D.  Camp 
Marshall  Space  Flight  Center,  Alabama  35812 

Dr.  F.  de  Percin 
Chief,  Atmospheric  Sci  BiMch 
Environments  Sciences  Division 
US  Army  Research  Office 
3046  Columbia  Pike 
Arlington,  VA  22204 

Geophysics  Officer 
Code  3260 

Pacific  Missile  Range 
Point  Mugu,  CA  93042 

US  Army  Liaison  Of^ce 
MIT,  Bldg  25,  Rm  131 
77  Massa^usetts  Avenue 
Cambridge,  Mass  02139 

The  Library  of  Congress 
ATTN:  Exchange  Sc  Gift  Div 
Washington,  D.  C.  20540 
2 

Head,  Atmospheric  Sci  Section 
National  Science  Foundation 
1600  G  Street  N.W. 

Washington,  D.  C.  20550 

Director,  S^tems  RAD  Service 
Federal  Aviation  Admin 
600  Indepexulence  Avenue  S.W. 

Washington,  D.  C.  20690 

Technical  Processes  Br  -  D628 
Room  806,  Libraries  Div  -  NOAA 
8060  Idth  Street 
Silver  Spring,  MD  20910 

William  A.  Main 
USDA  Forest  Service 
215  Natural  Resources  Building 
Michigan  State  University 
East  Lansing,  Mich  48823 

NASA  Headquarters 
Meteorology  &  Sounding  Branch 
(CODE  SAM) 

Space  Applications  Programs 
Washington,  D.  C.  20646 

Chief,  Fallout  Studies  Branch 
Division  of  Biolc^  Sc  Medicine 
Atomic  Energy  Commission 
Washington,  D.  C.  20546 

Director 

Atmospheric  Physics  &  Chem  Lab 
Room  31 

NOAA  Drat  of  Commerce 
B<Hilder,  Colorado  80302 

Natl  Center  for  Atmos  Research 
NCAR  Library,  Aoquisitions-Rept 
Boulder,  Colorado  80302 


E&R  Center,  Bureau  of  Reclamation 
ATTN:  D-751,  Bldg  67 
Denver,  Colorado  80225 

Natl  Oceanographic  Data  Center  -  D722 
Natl  Oceanic  &  Atmos  Admin 
US  Dept  of  Commerce 
Rockville,  MD  20852 

Institute  of  Sci  &  Technology 
The  University  of  Michigan 
Box  618,  (IRIA  Library) 

Ann  Arbor,  Mich  48107 

NASA  Scientific  &  Tech  Info  Fac 
ATTN:  Acquisitions  Br  (S-AK/DL) 
P.  O.  Box  33 
College  Park,  MD  20740 
2 

Battelle  Memorial  Institute 
Strategic  Technolo^  Office 
Information  Analysis  Center 
505  King  Avenue 
Columbus,  Ohio  43201 

Infrared  Info  &  Analysis  Center 
University  of  Michigan 
Institute  of  Science  &  Technology 
Box  618 

Ann  Arbor,  Mich  48107 

VELA  Seismic  Info  Center 
University  of  Michigan 
Box  618 

Ann  Arbor,  Mich  48107 

Commanding  General 
US  Army  Elect  Command 
ATTN:  AMSEL-RDLNF 
Ft  Monmouth,  NJ  07703 

Commanding  Officer 

US  Army  Engineer  Topographic  Labs 

Ft  Belvoir,  VA  22060 

CO,  USA  Foreign  Sci  &  Tech  Center 
ATTN:  AMXST-ISI 
220  7th  Street,  NE 
Charlottesville,  VA  22901 

Commanding  General 
US  Army  Elect  Command 
ATTN:  AMSEL-BL 
Ft  Monmouth,  NJ  07703 
3 


Commanding  General 
US  Army  Elect  Command 
ATTN:  AMSEL-BL-AP 
Ft  Monmouth,  NJ  07703 

Commanding  General 
US  Army  Elect  Command 
ATTN:  AMSEL-BL-TE 
Ft  Monmouth,  NJ  07703 

Commanding  General 
US  Army  Elect  Command 
ATTN:  AMSEL-BL-FMA 
Ft  Monmouth,  NJ  07703 


Environmental  Protection  Agency 
Division  of  Meteorology 

Research  Triangle  Park,  North  Carolina  27711 
Commandant 

USArmy  Field  Artillery  School 
ATTN:  Dept  of  Target  Acquisition 
Ft  Sill,  OK  73503 

Commandant 

US  Army  Field  Artillery  School 
ATTN :  Met.  Division 
Ft  Sill,  OK  73503 

Armament  Dev  &  Test  Center 
ADTC(DLOSL) 

Eglin  AFB,  FL  32542 

Gap^  D.  Atkinson,  LtC  USAF 
Chief,  Technical  Services  Div 
DCS/ Aerospace  Sciences 
ATTN:  AWS/DNTI 
Scott  AFB,  IL  62225 

Director 

National  Security  Agcy 
ATTN:  TDL 

Ft  George  G  Meade,  MD  20755 

Commander,  Naval  Ship  Sys  Command 
Tech  Library,  RM  3  S-08 
Washington,  D.  C.  20360 

Commander 

Naval  Weather  Ser  Command 
Washington  Navy  Yard  (Bldg  200) 

Washington,  D.  C.  20390 

Commanding  Officer 
Harry  Diamond  Laboratory 
ATTN:  Library 
Washington,  D.  C.  20438 

Rome  Air  Devel  Center 

ATTN:  Documents  Library  (TDLD) 

Griffiss  AFB,  NY  13440 

Air  Force  Cambridge  Rsch  Labs 
ATTN:  LC(1124) 

L.  G.  Hanscom  Field 
Bedford,  Massachusetts  01730 

Air  Force  Avionics  Laboratory 
ATTN:  AFAL/DOT,  STINFO 
Wright-Patterson  AFB,  Ohio  45433 
2 

HQ,  AFCS(EICRW) 

Richards-Gebaur  AFB,  MO  64030 

Ofc  of  Asst  Ch  of  Staff  for  Intelligence 
Dept  of  the  Army 
ATTN:  ACSI-DD 
Washington,  D.  C.  20310 
2 

Commanding  (xeneral 
US  Army  Materiel  Command 
ATTN:  AMCRD-H 
Washington,  D.  C.  20315 


Commanding  General 
US  Army  Missile  Command 
ATTN:  AMSMI-RR,  Dr.  J.  P.  Hallowes 
Redstone  Arsenal,  Alabama  35809 

CG,  US  Army  Missile  Cmd 
Redstone  Scientific  Info  Center 
ATTN:  Chief,  Document  Section 
Redstone  Arsenal,  Alabama  35809 
2 

Commanding  General 
US  Army  Test  &  Eval  Cmd 
ATTN:  NBC  Directorate 
Aberdeen  Proving  Ground,  MD  21005 

Commanding  Officer 
USACDC  FAA 
Ft  Sill,  OK  73503 

Commanding  Officer 
Frankford  Arsenal 
ATTN:  Library,  H1300,  Bl.  51-2 
Philadelphia,  PA  19137 

Conunanding  Officer 
Frankford  Arsenal 
ATTN:  WlOOO-65-1  (Mr.  Helfrich) 
Philadelphia,  PA  19137 

Commanding  Officer 
USASA  Test  &  Eval  Cen 
Ft  Huachuca,  Arizona  85613 

Commanding  Officer  &  Dir 
US  Navy  Ele  Lab 
ATTN :  Library 
San  Diego,  CA  92152 

Commander 

US  Naval  Ordnance  Lab 

ATTN :  Tech  Library 

White  Oak,  Silver  Spring,  MD  20910 

Commander 

Naval  Elect  Sys  Cmd  HQ 
ATTN:  CODE  05611 
Washington,  D.  C.  20360 

Commandant,  Marine  Corps 
HQ,  US  Marine  Corps 
ATTN:  CODE  A04C 
Washington,  D.  C.  20380 

Marine  Corps  Dev  &  Educ  Cmd 
Development  Center 
ATTN:  C-E  Div 
Quantico,  VA  22134 

Commander 

US  Naval  Weapons  Lab 
ATTN :  KEB-2F  (FENN) 

Dahlgren,  VA  22448 


Commander,  Naval  Air  Sys  Command 
Meteorological  Div  (AIH-540) 
Washington,  D.  C.  20360 

Air  Force  Cambridge  Rsch  Labs 
ATTN:  LKI 
L.  G.  Hanscom  Field 
Bedford,  Mass.  01730 
2 

AFCRL  (LYV) 

L.  G.  Hanscom  Field 
Bedford,  Mass.  01730 

HQ,  Elect  Sys  Div  (TRI) 

L.  G.  Hanscom  Field 
Bedford,  Mass.  01730 

Recon  Central/RSA 
AF  Avionics  Lab 

Wright-Patterson  AFB,  Ohio  45433 

HQ,  Air  Weather  Service 

ATTN:  AWVAS/TF  (R.  G.  Stone) 

Scott  AFB,  IL  62225 

US  Air  Force  Security  Ser 

ATTN:  TSG 

San  Antonio,  TX  78241 

HQ,  Air  Force  Sys  Command 
ATTN:  DLTE 
Andrews  AFB,  MD  20331 

Air  Force  Weapons  Lab 
ATTN:  WLIL 
Kirtland  AFB,  NM  87117 

Ofc  of  Asst  Ch  of  Staff  for  DS-SSS 
Dept  of  the  Army 
Rm  3C466,  The  Pentagon 
Washington,  D.  C.  20315 

Ofc  Asst  Sec  of  the  Army  (R&D) 
ATTN:  Asst  for  Research 
Rm  3-E-379,  The  Pentagon 
Washington,  D.  C.  20310 

Chief  of  Rsch  &  Development 
ATTN:  CRDES 
Washington,  D.  C.  20315 
2 


Commanding  General 
US  Army  Materiel  Command 
ATTN:  AMCMA-EE 
Washington,  D.  C.  20315 

Commanding  General 
US  Army  Materiel  Command 
ATTN:  AMCRD-TV 
Washington,  D.  C.  20315 
2 

Commanding  General 
US  Army  Missile  Command 
ATTN:  AMSMI-RRA,  Bldg  7770 
Redstone  Arsenal,  AL  35809 


CG,  US  Anny  Missile  Command 
Redstone  Scientific  Info  Cen 
ATTN:  Chief,  Document  Section 
Redstone  Arsenal,  AL  35809 
3 


USA  Security  Agcy  Combat  Dev  Actv 
ATTN:  lACDA.P(T)  &  lACDA-P(L) 
Arlington  Hall  Station,  Bldg  420 
Arlington,  Virginia  22212 
2 


Commanding  General 

US  Army  Combat  Devel  Command 

Intelligence  &  Control  Sys  Gr. 

ATTN:  INCSCSD-C 
Ft  Belvoir,  VA  22060 
2 

Commanding  Officer 
US  Army  Limited  War.  Lab 
ATTN:  CRDLWLTC 
Aberdeen  Proving  Ground,  MD  21005 

Director 

US  Army  Ballistic  Lab 

ATTN:  AMXRD-BSP  (Dr.  Reuyl) 

Aberdeen  Proving  Ground,  MD  21005 

Commanding  General 
US  Army  Test  &  Eval  Command 
ATTN:  AMSTE-EL 
Aberdeen  Proving  Ground,  MD  21(X)5 

Commanding  General 
US  Army  Test  &  Eval  Command 
ATTN:  AMSTE-FA 
Aberdeen  Proving  Ground,  MD  21(X)5 

Commanding  General 
US  Army  Test  &  Eval  Command 
ATTN:  AMSTE-NBC 
Aberdeen  Proving  Ground,  MD  21(X)5 

Commanding  General 
US  Army  Munitions  Cmd 
ATTN:  AMSMU-RE-R 
Dover,  NJ  07801 

Director 

US  Army  Munitions  Command 
Operations  Research  Group 
Aberdeen  Proving  Ground,  MD  21010 

Commanding  C^neral 
US  Army  Munitions  Command 
ATTN:  AMSMU-RE-P 
Dover,  NJ  07801 


Technical  Support  Directorate 
ATTN :  Technical  Library 
Bldg  3330 

Edgewood  Arsenal,  MD  26010 
Commandant 

US  Army  Chemical  Cen  &  School 
Micromet.  Section  (Chem.  Br.) 

Ft  McClellan,  Ala  36201 

Director 

USA  Engr  Waterways  Exper  Sta 
ATTN:  Ptsch  Center  Library 
Vicksburg,  Mississippi  39180 
2 

CG,  Deseret  Test  Center 
ATTN:  STEPD-TT.ME(S)  MET  DIV 
Bldg  103,  Soldiers  Circle 
Fort  Douglas,  Utah  84113 

Commanding  General 

USA  CDC  Combat  Arms  Groim 

Ft.  Leavenworth,  Kansas  66027 

Commanding  Officer 
Eustis  Directorate 
US  Army  Air  Mobility  R&D  Lab 
ATTN;  Technical  Library 
Fort  Eustis,  VA  23604 

Director,  US  Army  Advanced  MaU  Concepts 
Agcy 

ATTN:  AMXAM 
Washington,  D.  C.  20315 

Commanding  Officer 
USACDC  CBR  Agcy 
Ft.  McClellan,  Ala  36201 

Commanding  General 
US  Army  Materiel  Command 
ATTN:  AMCRD  R  (H.  Cohen) 

Washington,  D.  C.  20315 

Commanding  Officer 
US  Army  Combat  Dev  Cmd 
Communications-Elect  Agcy 
Ft  Monmouth,  NJ  07703 


Commanding  Officer 
Pica  tinny  Arsenal 
ATTN:  SMUPA-TVI 
Dover,  NJ  07801 

Commanding  Officer 
USA  Garrison 
ATTN;  Tech  Ref.  Div. 

Fort  Huachuca,  Arizona  85613 

Chief,  A.M.  &  EW  Division 
ATTN:  USAEPG-STEEP-TD 
Fort  Huachuca,  Arizona  85613 

US  Army  Rsch.  Off.  •  Durham 
Box  CM,  DUKE  Station 
ATTN:  CRDARD-IP 
Durham,  N.  Carolina  27706 


US  Army  Liaison  Office 
MIT-Lincoln  Lab,  Rm  A-210 
P.  O.  Box  73 
Lexington,  Mass.  02173 

HQ,  US  Army  Combat  Dev  Cmd 
ATTN:  CDCLN-EL 
Ft  Belvoir,  VA  22060 

Commanding  General 
US  Army  Tank- Automotive  Cmd 
ATTN :  A  MST A-Z,  M  r.  R.  McGregor 
Warren,  Michigan  48090 

Commanding  General 
US  Army  Elect  Command 
ATTN:  AMSEL-GG-TD 
Ft  Monmouth,  NJ  077Q3 
2 


Commanding  General 
US  Army  Elect  Command 
ATTN:  AMSEL-EW 
Ft  Monmouth*  NJ  07703 

Commanding  General 
US  Army  Elect  Command 
ATTN:  AMSEL-ME-NMP-PS 
Ft  Monmouth*  NJ  07703 

Commanding  General 
US  Army  Elect  Command 
ATTN:  AMSEL-MS-TI 
Ft  Monmouth*  NJ  07703 


Director*  Defense  Nuclear  Agency 
ATTN :  Tech  Library 
Washington*  D.  C.  20^5 

Naval  Ships  Systems  Command 
ATTN;  CODE  20626  (Tech  Lib) 
Main  Navy  Bldg.  Rm  1528 
Washington*  D.  C.  20325 

Director 

US  Naval  Rsch  Laboratory 
ATTN:  CODE  2627 
Washington*  D.  C.  20390 
2 


Commanding  General 
US  Army  Elect  (Command 
ATTN:  AMSEL-XL-D 
Ft  Monmouth,  NJ  07703 


CO,  Edgewood  Arsenal 

Defense  Devel  &  Engineering  Lab 

ATTN:  8MUEA-DDW-2,  Mr.  H.  Tannenbaum 

Edgewood  Arsenal*  MD  21010 


Commanding  General 
US  Army  Elect  Command 
ATTN:  AMSEL-WL-D 
Ft  Monmouth*  NJ  07703 


Commanding  Officer 
Aberdeen  ProviM  Ground 
ATTN:  STEAP-TL 
Aberdeen  Proving  Groimd*  MD  21002 


Commanding  General 
US  Army  Elect  Command 
ATTN:  AMSEL-NL-D 
Ft  Monmouth*  NJ  07703 


Director*  Ballistic  Rsch  Labs 
US  Array  Aberdeen  Rsch  &  Dev  Cen 
ATTN:  AMXRD-BTL*  Mr.  F.  J.  Allen 
Aberdeen  Proving  Ground,  MD  21005 


Commanding  Genet al 
US  Army  Elect  Command 
ATTN:  AMSEL-KL-D 
Ft  Monmouth*  NJ  07703 

Commanding  General 
US  Army  Elect  Command 
ATTN;  AMSEL-VL-D 
Ft  Monmouth,  NJ  07703 

Commanding  General 
US  Army  Elect  Command 
ATTN:  AMSEL-CT-D 
Ft  Monmouth*  NJ  07703 
3 


Conunanding  Officer 
US  Army  B^listic  Rsch  Labs 
ATTN:  AMXBR-B  &  AMXBR-IA 
Aberdeen  Proving  Ground*  MD  21005 
2 

Commanding  Officer 
US  Army  Land  Warfare  Lab 
ATTN:  CRDLWL-4A 
Aberdeen  Proving  Ground*  MD  21005 

Commanding  Officer 
US  Army  Elect  Proving  Ground 
ATTN:  STEEP-MT-EA 
Ft  Huachuca*  Arizona  85613 


Commanding  General 
US  Army  Elect  Command 
ATTN:  AMSEL-SC 
Ft  Monmouth*  NJ  07703 


d^ommanding  Officer 

US  Army  Dugway  Proving  Ground  Library 
ATTN:  STEDP-TL*  Tech  Library 
Dugway*  Utah  84022 


Commanding  General 
US  Army  Elect  Command 
ATTN:  AMSEL-RD-D 
Ft  Monmouth*  NJ  07703 

Commanding  General 
US  Army  Combat  Dev  Cmd 
ATTN:  CDCMR-E 
Ft  Belvoir,  VA  22060 


Commanding  Officer 
US  Army  Cold  Regions  R&E  Lab 
ATTN:  CRREL-RP  (Dr.  Yin  Chao  Yen) 
Hanover*  New  Hampshire  03755 

Commanding  Officer 
Yuma  Proving  Ground 
ATTN:  STEYP-AD(TECH  LIBRARY) 
Yuma*  Arizona  85364 


President 

US  Army  Artillery  Board 
Ft  Sill*  OK  73503 

Dir  of  Defense  Rsch  &  Engin 
ATTN :  Tech  Library 
RM  3E-1039*  The  Pentagon 
Washington*  D.  C.  20301 


Commanding  Officer 
US  Army  Arctic  Test  Center 
APO  Seattle  98733 
2 

(X),  US  Army  Tropic  Test  Center 
ATTN:  STETC-AD-TL 
APO  New  York  09827 
2 


Defense  Intelligence  Agcy  Commanding  Officer 

ATTN:  DIAAP-10A2  USACDC  Inst  of  Nuclear  Studies 

Washington*  D.  C.  20301  Ft  Bliss*  TX  79916 


Commamding  Officer 
USADDC  Inst  of  Systems  Analysis 
ATTN;  CDISA-SSCI  &  CDISA-STANO 
Ft  Belvoir,  VA  22060 

Commanding  General 
US  Army  Missile  Commamd 
ATTN:  AMSMI-RFG,  Mr.  N.  Bell 
Redstone  Arsenal,  AL  35809 

Commanding  Officer 
Frankford  Arsenal 

ATTN;  SMUFA.N-6300.  Mr.  M.  Schoenfield 
Bldg  201 

Philadelphia.  PA  19137 

Commanding  Officer 
US  Army  Satellite  Comm  Agcy 
ATTN;  AMCPM-SC-3 
Ft  Monmouth,  NJ  07703 

Director 

Night  Vision  Lab  (USAECOM) 

ATTN:  ASMEL-NV-OR,  Mr.  S.  Segal 
Ft  Belvoir,  VA  22060 

Sylvania  Elect  Sys-Western  Div 
ATTN:  Technical  Reports  Lib 
P.  O.  Box  205 
Mountain  View,  CA  94040 

NASA  HEADQUARTERS 
Prg.  Ch,  MET  &  Soundings  (SRM) 

Earth  Observations  Programs 
Washington,  D,  C.  20546 

Library-R51-Tech  Reports 
Environmental  Research  Labs 
NOAA 

Boulder,  Colorado,  80302 


Commamdamt 

US  Army  Air  Defense  School 
ATTN:  C&S  Dept,  MSL  SCI  DIV 
Ft  Bliss.  TX  79916 

Defense  Documentation  Center 
ATTN:  DDC-TCA 
Cauneron  Station  (Bldg  5) 

Alexamdriai,  VA  22314 
12 

Director,  USN  Research  Lab. 

ATTN:  Code  6530,  H.  Shenker 
Washington,  D.  C.  20390 

Mr.  E.  F.  Corwin 
Head,  Met  Research  Bramch 
Meteorological  Division 
Naval  Air  Sys  Cmd  (AIR-5401) 

Washington,  D.  C.  20360 

Commanding  Officer 
US  Army  Combat  Dev  Cmd 
Intelligency  Agency 
ATTN:  CDCINT-P,  Cpt  Thoresen 
Ft.  Huachuca,  Arizona 

Director 

US  Army  Advamced  Materiel  Concepts  Agency 
2461  Eisenhower  Avenue 
Alexandria,  VA  22314 

University  of  Okladioma 
Research  Institute 
Field  Artillery  Research  Off 
P.  O.  Box  3124 
Ft.  Sill,  OK  73503 

Dir  US  Naval  Research  Lab 
ATTN:  CODE  6530,  H.  Shenker 
Washington,  D.  C.  20390 


Air  Force  Weapons  Lab 
ATTN:  WLIL 
Kirtland  AFB,  NM  87117 
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